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We investigate Second Harmonic Generation (SHG) in monolayer WS 2 both deposited on a Si0 2 /Si 
substrate or suspended using transmission electron microscopy grids. We find unusually large second order 
nonlinear susceptibility, with an estimated value of d eff —4.5 nm/V nearly three orders of magnitude larger 
than other common nonlinear crystals. In order to quantitatively characterize the nonlinear susceptibility of 
two-dimensional (2D) materials, we have developed a formalism to model SHG based on the Green's 
function with a 2D nonlinear sheet source. In addition, polarized SHG is demonstrated as a useful method to 
probe the structural symmetry and crystal orientation of 2D materials. To understand the large second order 
nonlinear susceptibility of monolayer WS 2 , density functional theory based calculation is performed. Our 
analysis suggests the origin of the large nonlinear susceptibility in resonance enhancement and a large joint 
density of states, and yields an estimate of the nonlinear susceptibility value d eff = 0.77 nm/V for monolayer 
WS 2 , which shows good order-of-magnitude agreement with the experimental result. 



Two-dimensional (2D) materials, exhibiting extraordinary and novel properties not available in their bulk 
forms, are at the forefront of nanomaterial technologies 1 " 8 . For instance, the transition metal dichalcogenide 
(TMD) family of monolayer materials (such as MoS 2 , WS 2 ), unlike their bulk counterparts, possesses direct 
band gaps 9 " 12 and could thus generate strong photoluminescence 13 " 17 . The TMD family of monolayer 
materials have also garnered interest in nonlinear optics 1819 . Characterization of exfoliated WS 2 20 and synthesized 
MoS 2 monolayers 21 " 23 has revealed large second-order nonlinear susceptibility, as high as nm/V 22 ' 23 . Here we 
investigate SHG in monolayer WS 2 synthesized by chemical vapour deposition. We have developed a formalism 
to model SHG in 2D materials using the Green's function with a 2D nonlinear sheet source, which also takes into 
account the focused excitation geometry and the substrate effect as opposed to a simple plane-wave model. SHG 
in WS 2 monolayers both suspended and deposited on Si0 2 /Si substrates is studied. The estimated second order 
nonlinear susceptibility is approximately three orders of magnitude higher than that of common nonlinear 
optical crystals such as BBO 24 . Polarized SHG in WS 2 is studied to probe the tensorial nonlinear susceptibility, 
which is related to the structural symmetry of the 2D crystal. To gain further insight into the large nonlinear 
susceptibility, density functional theory (DFT) calculations of the second-order nonlinear susceptibility are also 
performed. 

Results 

A WS 2 monolayer lacks inversion symmetry and has non- vanishing second- order nonlinear susceptibility. An 
experimental setup was developed to characterize the SHG in the synthesized WS 2 samples (Fig. la). Fig. 1(b) 
shows an optical micrograph of a triangular monolayer of WS 2 grown on a Si0 2 /Si substrate. The same triangle is 
also shown as an SHG image in Fig. 1(c), obtained by raster scanning the sample with a step size of 1 um. As can 
be observed from the SHG image, the background contribution from the substrate surface in the surrounding area 
is negligible. To further isolate substrate effects, the same experiments were also performed on a WS 2 sample that 
was suspended on a transmission electron microscopy (TEM) grid. Fig. 1(d) shows an optical micrograph of the 
WS 2 film on the TEM grid, and Fig. 1 (e) shows the SHG image of the same triangle. Although the WS 2 triangle can 
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Figure 1 | (a) Schematic diagram of the experimental setup for 
characterizing SHG in WS 2 monolayers. The fundamental beam is from a 
mode-locked Ti:sapphire laser (832 nm center wavelength, 88 MHz 
repetition rate) and the second harmonic signal generated by the WS 2 
monolayer is epi-collected and measured using a spectrometer, (b) Optical 
micrograph of a synthesized triangular WS 2 monolayer grown on a Si0 2 /Si 
substrate, (c) SHG image of the same WS 2 triangle obtained by raster 
scanning the sample, (d) Optical micrograph of a synthesized triangular 
WS 2 monolayer placed on a transmission electron microscopy grid, (e) 
SHG image of the same WS 2 triangle obtained by raster scanning the 
sample. 
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Figure 2 | (a) SHG scan with sample translated along the axial direction. 
The distance is relative to the surface of the substrate. The vertical axis 
shows the SH signals generated from WS 2 on a Si0 2 /Si substrate (average 
incident power = 1 mW) and suspended on a TEM grid (average incident 
power = 0.7 mW), as well as SH signals for the bare Si0 2 /Si substrate and 
TEM grid, by focusing the fundamental beam to different depths. There is 
no detectable SHG signal from the bare substrate without the WS 2 
triangular monolayer, thus allowing all the SH signal generated at the 
surface of the sample to be attributed to the presence of WS 2 . (b) Power 
dependence of the SH signal generated from a triangular WS 2 monolayer in 
a logarithmic scale on both the suspended and the Si0 2 /Si substrate. The 
measured power is shown as points and the linear fit is plotted as a solid 
line. The bottom inset shows the normalized fundamental laser spectrum 
and a typical normalized spectrum of the SH signal generated from a WS 2 
monolayer. 



barely be seen on the TEM grid micrograph due to weak reflection 
contrast, the triangle is easily apparent in the SHG image. A depth 
scan was also performed on the WS 2 samples and bare substrates to 
measure the SHG signal as a function of sample depth position, by 
translating the sample axially with a computer- controlled linear 
stage. The results shown in Fig. 2(a) indicate that the Second 
Harmonic (SH) signals generated from both the on-substrate and 
suspended samples are predominantly from the WS 2 monolayers, as 
the bare substrates had no detectable SH signal under the same 
measurement conditions. To characterize the second-order non- 
linear susceptibility, we measured the average power of the SHG 
signal as the power of the fundamental beam was varied. Fig. 2(b) 
shows that the resulting relation is quadratic: a linear fit of the log-log 
plot reveals a slope of 2.04 for the TEM-suspended WS 2 sample, and 
2.12 for the WS 2 /Si0 2 /Si substrate sample. The lower inset of 
Fig. 2(b) shows representative spectra of the fundamental and SHG 
signal, which also clearly indicate the frequency doubling. 

Discussion 

In order to quantitatively characterize the second- order nonlinear 
susceptibility, we have developed a formalism to model SHG from 
2D materials on a semi-infinite substrate. A plane-wave fundamental 



beam uniformly filling the aperture of a lens (focal length F, aperture 
radius a, numerical aperture NA = a/F) is focused onto the 2D 
material, and the generated SH signal is epi-collected. The fun- 
damental field at the sample can be represented by the Debye integ- 
2 
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where E 10 represents the incident field at the lens aperture, n x is 
the refractive index of the substrate at the fundamental wavelength 
Xi, and circQ is the circular aperture function. The paraxial approxi- 
mation is used, and the Fresnel coefficient for normal incidence is 
also taken into account. The SHG is governed by the wave equation 
(S/ 2 + w 2 2 c~ 2 n 2 (z))E 2 = — n 0 a}\£ 0 e 2 ")^'. eieiE 2 (x,y,0), where E 2 is 
the complex SH field, e\ and e 2 represent the polarization of the 
fundamental and SH beam, co 2 is SH angular frequency, c is the speed 
of light in vacuum, e 0 and /n 0 are the vacuum permittivity and per- 
meability, and x (2) is the second-order nonlinear susceptibility ten- 
sor 19 . Due to its atomic thickness, the nonlinear susceptibility of the 
2D material can be modelled by a Dirac delta function with an 
effective surface nonlinear susceptibility, i.e. e 2 ~x^ '■ = xf^( z )- 
Additionally, n(z) = 1 (z < 0) or n 2 (z > 0), the refractive index of the 
substrate at the SH frequency. The SHG from the substrate surface is 
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neglected, as justified by the experimental observation of no signal 
for the substrate background measurements of Fig. 1(c) and Fig. 2(a). 

The relationship between the fundamental excitation power P x 
and the SH signal power P 2 is obtained by solving for the Green's 
function with the nonlinear sheet source: 



P 2 =- 



32NA 2 (xl 2) ) 2 Pi 2 
£ 0 c^(n 2 + l) 2 ( ni + l) 
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with 0 = 87T cos p 
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From Eq. (1) we can express the surface nonlinear susceptibility in 
terms of the average power, the pulse width, and the pulse repetition 
rate: 



2 Jaf(H 2 + l) 2 (»i + l) 4 



32NAH 2 P 2 avl <l> 



(2) 



where R is the repetition rate, t { is the pulse width, and P avi is the 
average power (i = 1: fundamental, 2: SH). Since the WS 2 monolayer 
has atomic thickness much shorter than the pulse length involved, 
the continuous wave approximation is still valid even though chirped 
femtosecond pulses are used in our experiment. In our calculation, 
the instantaneous power of the fundamental beam is approximated 
by the pulse energy divided by the pulse duration time. The calcu- 
lated nonlinear susceptibility represents an averaged value over the 
pulse bandwidth. Based on data taken similar to that of Fig. 2(b) 
(from 30 triangles on 3 separate substrates) and Eq. (2), the effective 
bulk-like second-order susceptibility of the grown WS 2 monolayers, 
with d e ff = xf)/2T where T = 0.65 nm is the thickness of a WS 2 
monolayer, was estimated to be 4.46 nm/V for the suspended film 
with a standard deviation of 0.09 nm/V across all samples, and 
4.51 nm/V for the WS 2 triangle on a Si0 2 /Si substrate with a stand- 
ard deviation of 0.55 nm/V. Note that the standard deviations should 
be interpreted as the consistency of the measured SHG power rather 
than the accuracy of the susceptibility value due to the assumptions 
that were made in our calculation. These estimated values are nearly 
three orders of magnitude greater than that of typical nonlinear 
crystals 24 . It should also be noted that the Si0 2 /Si substrate may 
not be simply treated as a semi-infinite homogeneous medium due 
to the presence of the Si0 2 /Si interface, as it is well known such 
substrate has contrast enhancement 26 . The calculated nonlinear sus- 
ceptibility thus needs to be adjusted by an enhancement factor, i.e., 
xf\ =tfmi£s ) where and r| 2 are the field enhancement factors 
for the fundamental and the epi- detected SH field, respectively, due 
to the Fabry-Perot cavity formed by air/Si0 2 and Si0 2 /Si interfaces. 
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Figure 3 | Polarized SHG in WS 2 triangles. The SHG images show two 
monolayer WS 2 triangles on the same Si0 2 /Si substrate. The polar plots 
show the normalized polarization dependence of SHG in WS 2 monolayer 
islands. Dots are measured data and the solid curves are the theoretical fit. 
The two polar plots reveal a rotation of 18 degrees between the two 
triangles. 



is the calculated value using Eq. (2), and x^ AT is the actual non- 
linear susceptibility. For a substrate with 300 nm of Si0 2 above sil- 
icon, the susceptibility value is estimated to be 15% more when the 
cavity enhancement factors are accounted. This analysis of the 
enhancement factor assumes a normally incident plane wave. A 
more rigorous approach would require angular spectrum decom- 
position of the incident wave. For this reason, the suspended sample 
is likely to yield a better estimate of the nonlinear susceptibility. 

In addition to the value of d eff , the tensorial property of the second- 
order nonlinear susceptibility of WS 2 monolayers can also be probed 
by polarized SHG. The WS 2 monolayer has P6 3 /mmc crystal sym- 
metry, resulting in a second- order nonlinear susceptibility tensor 
with nonzero elements dyyy = — d yxx = — d xxy = — d xyx where x 
and y represent the crystal axes 19 . The polarized SH signal power 
P 2 oc cos 2 30 where 0 is the angle between the x crystal axis and 
the incidence polarization 27 . This relationship is shown compared to 
experimental results in Fig. (3). Two WS 2 monolayer triangles on the 
same Si0 2 /Si substrate are shown by SHG raster scan. The two inserts 
are polar plots where the circles are measured data points and the 
solid line is the cos 2 30 fit. The two polar plots are rotated by 18 
degrees with respect to each other, and the left edge of the two 
triangle SHG images appear to be rotated by a similar angle. The 
experimental data and theoretical curve agree well and clearly dem- 
onstrate six-fold symmetry. As the tensor structure of x (2) is closely 
related to the symmetry of the material structure, polarized SHG can 
be a useful tool for probing the structural symmetry of 2D materials. 

To understand the origin of the giant experimentally observed 
second-order nonlinearity we have also calculated the second-order 
nonlinear optical susceptibility of monolayer WS 2 at the level of 
density functional theory (DFT), extending beyond existing literat- 
ure results in the tight binding approximation 28 . Note that the band 
gap measured by photoluminescence agrees with the DFT calcula- 
tion surprisingly well, due to the partial cancellation of quasi -particle 
effects and exciton binding 29 . On the other hand, calculations at the 
quasi-particle or Bethe-Salpeter equation level are computationally 
expensive, due to the very fine k-point mesh required. For an accur- 
ate description of the second-order response we follow the formalism 
of Sipe and Ghahramani 30 33 as implemented in the ABINIT pack- 
age 34 ' 35 , including both inter-band transitions and intra-band cur- 
rents. Although completely filled bands produce no intra-band 
current at the linear order, intra-band motion does occur at second 
order 30,32 . The ground-state properties and response functions were 
calculated within the local density approximation 35 using norm- 
conserving pseudopotentials with an energy cut-off of 760 eV and 




Figure 4 | The band structure of monolayer WS 2 . The parallel bands along 
T-K give rise to maxima in the joint density of states. 
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Figure 5 | The nonlinear susceptibility/^ of monolayer WS 2 . Also shown 
is the colour-coded joint density states contributions from the highest 
valence band to the first, second, and third conduction bands, aligned to 
the calculated spectra at the appropriate second-order resonant 
frequencies and following the colour scheme of Fig. 4. 

a k-point mesh of 69 X 69 X 1 within ABINIT 36 " 40 . The optical 
spectra achieved convergence with 50 conduction bands and a 
smearing of 15 meV. 

The band structure of monolayer WS 2 features three nearly par- 
allel bands along T-K, as shown by the shaded areas near II and IV in 
Fig. 4. The highest valence band is parallel to both the first and the 
second conduction bands, yielding two peaks labelled II and IV in 
the joint density of states (JDOS) of Fig. 5. Another prominent fea- 
ture in the JDOS, labelled V, originates from transitions near the K 
point. The second- order interband transition V could be strongly 
enhanced by the intermediate transition I at the half-way point in 
energy, if these matrix elements are non -vanishing. This is confirmed 
by the calculated nonlinear spectra. 

In addition to the calculated # (2) , Fig. 5 also shows the JDOS, 
horizontally rescaled to j(2(o) to match resonances arising from the 
2o) terms and separated out by colour-coded destination bands fol- 
lowing the same colour scheme as in Fig. 4. The main features in the 
JDOS match well with the calculated x (2) spectra. Resonance I has the 
lowest excitation energy, corresponding to half the band gap. 
Resonances II and IV clearly benefit from a large JDOS. Resonance 
V arises from transitions near the K point, which lies at almost twice 
the band gap (also at K). Near transition II (corresponding to a 
wavelength of 832 nm) we obtain / (2) = 0.5 nm/V. With the simu- 
lation cell height of 2 nm, we obtain a value of d eff =0.77 nm/V. The 
calculation confirms the very large non-linear response observed 
experimentally, and considering the overall magnitude of the effect 
and the neglect of quasiparticle and excitonic effects, the semi-quant- 
itative agreement to experiment is surprisingly good. 

In summary, we report strong SHG from synthesized WS 2 mono- 
layer islands both on Si0 2 /Si substrate and suspended on a TEM grid, 
as well as a theoretical calculation of the second- order susceptibility 
using DFT. In order to experimentally determine the second-order 
susceptibility of 2D materials, we have developed a SHG model based 
on the Green's function, yielding an estimated bulk-like second- 
order nonlinear susceptibility d eff approximately 4.5 nm/V for a 
WS 2 monolayer. In addition, by using polarized SHG, the crystal 
symmetry and orientation of the WS 2 monolayer can be revealed, 
demonstrating SHG as a useful method to probe structural informa- 
tion in 2D materials. Further analysis through DFT calculations 
indicate that the large second- order susceptibility is due to resonance 
enhancement as well as the large joint density of states. Despite the 
approximations involved, the calculated nonlinear susceptibility (d eff 
= 0.77 nm/V) using DFT shows reasonable order- of- magnitude 
agreement with the experimental result. With giant / (2) and sub- 
nanometre thickness, the integration of 2D TMD materials with 
photonic circuits to realize new nonlinear optical devices could be 
a fertile ground worth further exploration. 



Methods 

Synthesis. The WS 2 triangle monolayer samples were grown on Si0 2 /Si substrates 
using a two-step process 17 . Briefly, one-nm films of W0 3 (99.998%, Alfa Aesar) were 
thermally evaporated onto a Si0 2 /Si substrate at 10" 6 Torr. The films were placed 
into a quartz reaction tube next to a boat containing sulphur powder (99.5%, Alfa 
Aesar). The sulphur zone was heated up to 250°C and the furnace where the W0 3 
sample was located was heated to 800°C. The synthesis was carried out in an inert 
environment under an Ar flow using 100 seem and atmospheric pressure at the outlet. 
To verify the structural quality and integrity of the grown WS 2 monolayers, Raman 
spectroscopy using a Renishaw in Via micro-Raman was carried out at 514.5 nm laser 
excitation. The 2LA(M) phonon mode (352 cm -1 ) had an intensity twice as high as 
that of the Alg mode (418 cm" 1 ), which, according to reference 41, is a characteristic 
of monolayer WS 2 , thus confirming that the synthesized WS 2 samples are indeed 
monolayers. 

Characterization. The output from a mode-locked Ti:sapphire laser (from KM Labs 
with centre wavelength of 832 nm and repetition rate of 88 MHz) was filtered, 
attenuated, and focused onto a sample by a long- working- distance objective lens 
(Mitutoyo 50 X, NA = 0.55, spot size —1.8 urn at fundamental wavelength). The 
generated SH signal was then back collected by the same lens, separated using a 
dichroic mirror, and filtered by a 405 nm bandpass filter before entering a 
spectrometer (PI Acton 2500i with a liquid nitrogen cooled charge coupled device — 
CCD camera). The signal value of the spectrometer was calibrated using an 
attenuated (by neutral density filters) SH signal separately generated from a BBO 
crystal with known power measured by a power meter (Newport 1830-C). In order to 
assist the alignment and ensure that the signal was collected from a WS 2 monolayer 
island, an imaging system was co-set up so that a removable mirror could direct the 
image signal towards an imaging CCD camera. In addition, in order to calculate the 
susceptibility, the laser pulse width needs be known in situ, i.e. directly on the sample 
at the focal point of the objective lens. Therefore, we performed a collinear frequency 
resolved optical gating (cFROG) measurement 42 by using the same WS 2 sample to 
determine the pulse width (106 fs, chirped). Although there exists dispersion in the 
X (2) of WS 2 (c.f. fig. 5), with our laser's relatively narrow bandwidth of about 18 nm, 
this was determined not to be an issue and was confirmed by the agreement of the 
retrieved spectrum with the measured laser spectrum. WS 2 monolayer may also 
enable the characterization of complex pulses with broad bandwidth if its 
susceptibility is pre-calibrated. To measure the polarization dependence, the 
fundamental beam was first prepared in a circular polarization state using a quarter- 
wave plate. A broadband polarizer was inserted between the objective lens and the 
dichroic mirror to serve as both the polarizer for projecting the input fundamental 
beam into a linear polarization state and the analyzer for the generated SHG signal. 
The polarizer was mounted on a computer-controlled rotational stage to 
continuously rotate its polarization axis. 
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